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Abstract — This paper investigates amplify-and-forward (AF) 
schemes for both one and two-way relay channels. Unlike most 
existing works assuming independent noise at the relays, we 
consider a more general scenario with correlated relay noise. 
We first propose an approach to efficiently solve a class of 
quadratically constrained fractional problems via second-order 
cone programming (SOCP). Then it is shown that the AF relay 
optimization problems studied in this paper can be incorporated 
into such quadratically constrained fractional problems. As a 
consequence, the proposed approach can be used as a unified 
framework to solve the optimal AF rate for the one-way relay 
channel and the optimal AF rate region for the two-way relay 
channel under both sum and individual relay power constraints. 

In particular, for one-way relay channel under individual relay 
power constraints, we propose two suboptimal AF schemes in 
closed-form. It is shown that they are approximately optimal in 
certain conditions of interest. Furthermore, we find an interesting 
result that, on average, noise correlation is beneficial no matter 
the relays know the noise covariance matrix or not for such 
scenario. Overall, the obtained results recover and generalize 
several existing results for the uncorrelated counterpart. 

I. Introduction 

Since the ampUfy-and-forward (AF) relay scheme was intro- 
duced, it has been widely studied in the context of cooperative 
communication ||T]-|[3l. It is an interesting technique from 
the practical standpoint because the complexity and cost of 
relaying, always an issue in designing cooperative networks, 
are minimal for AF relay networks. As the simplest coding 
scheme, AF is also used to estimate the relay network capacity. 
Obviously, the achievable rate of the AF scheme can be viewed 
as a lower bound to the network capacity. In addition to its 
simplicity, AF is known to be the optimal relay strategy in 
many interesting cases, e.g., ID and IS). 

In certain resource constrained networks, collaborative AF 
schemes have been developed to exploit the multi-antenna gain 
in multiple relay networks. The problem of finding the optimal 
AF scheme for two-hop one-way relay channel under sum 
or/and individual relay power constraints has been extensively 
studied in previous works ||6l-||9]. Under the sum power con- 
straint, an algebraic approach was proposed to find the optimal 
AF scheme in a closed form in terms of maximizing the 
transmission rate ID, ||9l. However, the problem becomes more 
challenging when the relays are subject to individual power 
constraints. For this problem, a semi-closed form solution has 
been derived in ||6]. The authors in Q also developed an 
iterative algorithm to solve the problem. In llTTI . the algebraic 



approach was generalized to solve the sum/individual rate 
maximization problem for a two-hop multiple access relay 
channel. Moreover, the AF scheme design issues have been 
considered for more general layered llT3l - [lT5l and non-layered 
[fT2| relay networks. 

Recently, analog network coding (ANC) ifTSl extends the 
AF-based one-way relaying scheme to a two-way relay chan- 
nel to support communications in two directions via a two- 
step protocol. The problem of finding the optimal rate region 
of a two-way relay channel with a two-step AF protocol has 
been studied for different network setups. The authors in llTTl 
characterized the maximum achievable rate region for the two- 
way relay channel with a single multi-antenna relay; while the 
maximum rate region for the AF two-way relay channel with 
multiple single-antenna relays has been studied in [18] and 
IIT9I under both sum and individual relay power constraints. 

The AF scheme designs mentioned above has been built 
upon the assumption of independent relay noise. However, in 
wireless networks, noise correlation between nodes may occur 
due to the common interference or noise propagation, and the 
previous design approaches cannot be directly applied to such 
correlated scenarios. For a two-hop relay network with noise 
correlation, it is worthwhile to mention the pioneer work llTol . 
which presented a closed-form solution under the sum relay 
power constraint. With the individual relay power constraints, 
although the algorithm in fl] can be generalized to solve this 
problem, its convergence rate is not guaranteed theoretically. 

In this paper, considering both sum and individual power 
constraints at the relays, our focus is on the design of optimal 
AF schemes for both one and two-way relay channels with 
correlated relay noise. Naturally, our design problem is more 
general and difficult than the existing ones. For simplicity, we 
assume that all the network parameters are real. Our main 
contributions are summarized as follows. 

1) We study a class of quadratically constrained fractional 
problems which is generally non-convex. With the aid 
of some transformation tricks, we show that this class of 
problems can be recast as a collection of convex SOCP 
problems and thus can be efficiently solved via interior 
point methods. The reason why we study such problems 
is that they can be used as a unified framework for 
formulating the AF relay optimization problems for both 
one and two-way relay channels under different power 
constraints. 
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Fig. 1. One-way Relay Channel 



2) We obtain the optimal AF scheme for one-way relay 
channel using the proposed formulation. Specifically, it 
is sufficient to get the optimal solution via solving only 
one SOCP problem. For individual power constraint, 
we also propose two suboptimal AF schemes which 
are approximately optimal under certain conditions of 
interest. Finally, we show that the noise correlation is 
beneficial to the AF performance on average. 

3) We also study the two-way AF relay channel without the 
assumption of channel reciprocality. We apply the pro- 
posed formulation rather than the non-convex weighted 
sum rate maximization (WSRMax) problem to find the 
boundary point of the union rate region. In particular, 
each boundary point can be calculated by solving two 
SOCP problems. 

Notation: diag {x} denotes a diagonal matrix with the ele- 
ments in X on its main diagonal and vec {X} = [xi, - ■ ■ , x„] 
with Xk on the diagonal of X. We use 7\A(/x,K) to denote 
n-dimensional joint Gaussian distribution with means and 
covariance matrix K. log(-) denotes the logarithm in the base 
2 and E[ ] denotes the expectation of a random variable. 

II. A Unified Framework 

As is shown in great many works, designing the optimal AF 
schemes in relay networks is formulated as an optimization 
problem. In this section, we provide a class of quadratically 
constrained fractional problems which can be used as a unified 
framework for formulating such problems for both one and 
two-way AF relay channels. 

Let's consider the quadratically constrained fractional prob- 
lem shown as follows. 

x^bobjfx 



max 



s.t. 



x^Aqx + Co 
x^A.x - x^b,b?'x - 



<0, 



X 



AjX + Cj < 0, j = mi + 1, ■ 



1, • • ■ mi ' 



(PI) 



where > 0, i = 0, 1, • • ■ mi, and A.; for i = 0, • • • , mi and 
A j for j = mi + l,-- - ,mi+m2are symmetric positive semi- 
definite matrices. It is observed that the objective function of 
the above problem is not concave and the first mi constraints 
are not convex. As a result, it is not a standard convex 
optimization problem. How to solve this problem is the focus 
of this section. With the aid of some transformation tricks, we 
recast (PI) as follows. 

Proposition 1: Problem (PI) can be converted into an 
equivalent quadratically constrained quadratic programming 



(QCQP) problem (P2) given as follows. 

max boy 

y,i' 

s.t. y^Aoy + c^v^ - 1 < 

y^Aiy - y^b.bf y + qu^ < 0, i = 1, • • • mi 

y'^ AjY + Cjv'^ < 0, j = mi + 1, • ■ • , mi + m2 

(P2) 

Proof: The proof is given in Appendix A. ■ 
It is first noticed that the objective function of (P2) is 
linear and thus concave. Then to solve (P2), we check the 
convexity of the feasible set. It is shown that except the ones 
in the second form, i.e., y-^A^y — y-^b^bf y + av'^ < for 
i = 1, • • • , mi, all the others are convex constraints. However, 
it is found that each constraint in the second form is a union 
of two convex second-order cone constraints. Consequently, 
the optimal solution of (P2) can be obtained by solving 
a collection of convex subproblems. It is concluded in the 
following proposition. 

Proposition 2: Problem (P2) can be efficiently solved by 
solving a collection of 2"'^ SOCP problems. 

Proof: The proof is given in Appendix B. ■ 
It should be pointed out that given the number mi, the 
computational complexity of the approach which solves (PI) 
is a polynomial function in the dimension of x. 

III. One-way Relay Channel 

In this section, we consider a two-hop one-way AF relay 
channel with a single source-destination pair {S, D) and n 
relay nodes as depicted in Fig. 1. It can be regarded as a 
layered relay network with three layers. Due to the channel 
fading, the direct path between the source and destination 
nodes is neglected. The channel coefficients are all real-valued 
numbers and remain constant during the operation. Therefore, 
we only focus on the design of the amplification gains at the 
relays in different power constraint setups. The full channel 
state information is revealed to all the network nodes. 

Assume that all the relay nodes are equipped with a single 
antenna and work in a half-duplex mode. Therefore, the data 
transmission takes place in two steps. During the first step, the 
source node sends xs with fixed transmit power Ps- The fcth 
relay receives yk = hs^kXs + Zk with Zk being the Gaussian 
noise and hs^k being the channel coefficient between source 
node S and relay node k. During the second step, the fcth 
relay sends Xk = akUk with ak being the amplification gain 
and hence the received signal at the destination node can be 
expressed as 



= a^^H2hiX5 

equivalent source signal 



Az 



ZD 
..T 



a H2Z + zjj 



(1) 



equivalent Gaussian noise 

,a„}, a ~ vec- 

= [hi.D, ■ • • , h,i,D\ 
channel vectors from the source node to the relays and 
from the relays to the destination node respectively, = 



where A = diagjai,- 
[hs,i, ■ ■ ■ ,hs,nf and h2 



jjA}, hi = 
^ denote the 



diag{hi}, for i ~ 1,2, and z = [zi,--- , Zn] and zjj are 
the Gaussian noise at the relays and the destination node 
respectively. The Gaussian noise is independent of the source 
signal. As the focus of the paper, we consider the Gaussian 
noise at the relays being correlated which is drawn according 
to 7V(0, K), but it is independent of zd ^ (O, a^). 

From ([T]i, it follows that the two-hop AF relay channel 
can be regarded as an equivalent point-to-point Gaussian 
channel where h^AhiXs is regarded as the source signal and 
h^Az + ZD is viewed as the Gaussian noise. As is well- 
known, the source node adopts the Gaussian codebook with 
xs drawn according to TV (0, Ps). The achievable rate is given 
by i? = 0.25 log (1 + SNR), where the pre-log factor is due 
to the half-duplex assumption and SNR denotes the received 
signal-to-noise ratio (SNR) at the destination node. Note that 
log( ) is an increasing function. Therefore, to maximize the 
transmission rate is equivalent to maximizing the received 
SNR at the destination node. 

We first consider a general case when the relays are 
subject to both sum and individual power constraints, i.e., 
{PsHj + K I) a < P,,„„ and al {a'l + hlj^Ps) < 
Pk for k = 1, • • • ,n, where Pgum is the sum power budget, 
Pk is the individual power budget for relay node k, is 
the kth diagonal element of K, and denotes the point- 
wise multiplication of the two matrix. Then the optimal AF 
scheme can be obtained via solving the following optimization 
problem in the form (PI). 

a^HahihfHaaPs 



max SNR (a) = 



s.t 



a^HaKHaO! + 

- K I) a < P,, 



(2) 



a. 



IfcCc < Pk, k = 1, - ■ ■ ,n 



where Ifc = diag {0, ■ • • , 0, 1, 0, • • ■ ,0} is a diagonal matrix 
with the fcth element equal to 1. 

By Proposition 1, it can be recast as 



max hi Ha^V^'s 

s.t. /3^H2KH2/3 + a^v^ - 1 < 

{PsUl + K I) /3 - P^urnv^ < 0- (3) 

fc = 1, • ■ • ,n 

Then by Proposition 2, to solve Q, only one SOCP problem 
should be solved, which is given as follows. 

max hfHs/Sy^ 

I3,v 



S.t. H2KH2/3 + a^v^ - 1 < 



(PsH2 + K0I)^/3||2 < /f^ 



(4) 



Psh^ 



S.k ^ 
1, • ■ • ,n 



Denote the optimal solution of ^ by (/3opt, Wopt), then the 
optimal solution of (|2]i is given by ccopt ~ /3opt/vopt- 



In particular, if K is a diagonal matrix, the above problem 
reduces to the case without relay noise correlation. Thus the 
proposed approach is also applied to the case. Although the 
closed-form solution cannot be obtained, we find the proposed 
approach efficient enough for application in practice. 

From the practical perspective, we are especially interested 
in the case when all the relays are only subject to the individual 
power constraint. This model can be used to characterize 
a network with each node having individual power supply. 
Moreover, to the best of our knowledge, the proposed approach 
is the most efficient one to solve the problem of finding the 
optimal AF scheme for the relay network with correlated relay 
noise and subject to individual relay power constraint. 

To further reduce the computational complexity, two subop- 
timal AF schemes are proposed and the application conditions 
of them are discussed in the sequel. Given the network 
parameters, it is easy to compute the following two schemes: 

1) Suboptimal scheme 1 (for relatively small cr^): 

ai =co(H2KH2)"'H2hi, (5) 
where ao = [ao.i, ■ • ■ ,ao,«]^ = (H2KH2)~^ H2hi, 



^, for k = 1, • • • , n, and co 



-,fc= 1,. 



2) Suboptimal scheme 2 (for relatively large a ): 



a.2 = [sign(/ii)ai^„ 



,sign(/i„)Q;„^„ 



(6) 



where hk = hs.khk.D^ for k = I,-- - ,n, and 
sign [hk) = 1 for /ife > otherwise sign [hk) = —1. 
Then we give the following definitions according to a.i, 
i = 1,2. 

Definition 1: We say that the relay noise dominates the sum 
noise at the destination node if 



Q:fH2KH2Q:i > 



(7) 



and that the destination noise dominates the sum noise at the 
destination node if 



a^H2KH2a2 < o-^, 



(8) 



Based on such definitions, we have the following results. 

Proposition 3: The achievable rate of oti is at most 0.25 
bit away from the corresponding optimal one when the relay 
noise dominates the sum noise, and the achievable rate of 0.2 
is also at most 0.25 bit away from the corresponding optimal 
one when the destination noise dominates the sum noise. 

Proof: The proof is given in Appendix C. ■ 

Finally, we verify our analytical results via the following 
example. 

Example: Let's consider a relay channel with three relay 
nodes. The network setups are given below: 



Ps 



1, = 1, hi = [l,2,lf , h2 = [3,3,1] 



K = 



1.5004 1.3293 0.8439 
1.3293 1.2436 0.6936 
0.8439 0.6936 1.2935 



TABLE I 

Optimal and suboptimal AF schemes 



Pk Scheme ai 02 03 R (ct) 

optimal 0.2000 0.1381 0.2088 0.2421 

suboptimali -0.1247 0.1318 0.0394 0.0779 
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suboptimal2 


0.2000 


0.1381 


0.2088 


0.2421 


optimal 


-0.8455 


0.9765 


0.4384 


0.8006 


suboptimali 


-0.8819 


0.9765 


0.2788 


0.7855 


suboptimal2 


1.4141 


0.9765 


1.4765 


0.3263 



Here we fixed the transmit power of the source node. For 
brevity, assume that all the relay nodes have the same individ- 
ual power budget. The optimal and suboptimal AF schemes 
are shown in table I, where R (a) = 0.25 log (1 + SNR (a)). 
From table I, we observe that when = 0.1, the correspond- 
ing achievable rate of 0.-2 is optimal and when Pk — 5, the 
corresponding achievable rate of a.i is approximately optimal. 
As a consequence, the numerical results coincide with the 
analysis given above. 

So far, we have studied the performance of the optimal and 
suboptimal AF schemes when the relay noise is correlated. We 
still curious about the effect of the correlation compared with 
the independent scenario. Does correlation help? This question 
has been answered in fTOl under sum relay power constraint 
scenario. However, it seems hard to answer the question under 
individual relay power constraint scenario due to the lack of 
analytical optimal AF scheme. Fortunately, it can be handled 
with the aid of our previous work in IITTI . To investigate the 
issue, we assume that all the channel coefficients are i.i.d 
Gaussian random variables with zero mean and unit variance. 
Then we obtain the following results. 

Proposition 4: For any source power P5 and relay power 
constraint Pk, k ~ I,-- - ,n, the AF performance under 
correlated relay noise is better than or equal to that under 
independent relay noise in terms of both average SNR and 
average rate over all the channel realizations. 

Proof: The proof is given in Appendix D. ■ 

IV. Two-way Relay Channel 

In this section, we consider a two-hop two-way AF relay 
channel consisting of two source nodes 5*1 and S2 and n relay 
nodes as shown in Fig. 2. We still assume that no direct paths 
between Si and S2 exists. The uplink channels from 6*1 and 
S2 to relay node k are denoted as /i ^ and f2.k respectively, 
while (71 fc and g2,k denote the downlink channels from relay 
node k to Si and S2 respectively. All the involved channels 
are assumed to take real values and remain constant during the 
operation period. In addition, all the channel state information 
is revealed to all the network nodes. Again, we assume the 
relay nodes receive correlated Gaussian noise which is drawn 
according to A/" (0, K). 

To achieve higher spectrum efficiency, a two-step AF proto- 
col is adopted, which is referred to as analog network coding 
lfT6l . During the first step, both source nodes simultaneously 



transmit their signals xsi and with fixed power Psi 
and Ps2 respectively to the relay nodes. The kth relay node 
receives yk = fi,kXSi + h^xs^ + Zk with Zk being the 
Gaussian noise. In the second step, the received signal at 
relay node k is multiplied by the amplification gain ak and 
retransmitted. After the self-interference cancelation, the two 
source nodes receive: 

yo,i = gfAfaXs^ + gfAz^+zo^ , (9) 

equivalent source signal equivalent Gaussian noise 
y0,2 = g2^A^^ + g^Az^+ Zo,2 , (10) 
equivalent source signal equivalent Gaussian noise 

where z = [zi, • • • , denotes the noise vector at the 
relays received in the first step, A = diagjai,-- - ,«„}, 
zo,i is the Gaussian noise with zero mean and variance ctq ^ 
received at Si in the second step, ~ • • • ,/i,n]"'" and 
gi = [.9i.ir'' ,ffi,n]"'", for i ~ 1,2. Assume that the relay 
noise z is independent of the source signals and the noise z^^i, 
i ~ 1,2. Similar to the one-way case, given the amplification 
gains, the two-way relay channel can be regarded as two point- 
to-point Gaussian channels. Therefore, both the source nodes 
adopt Gaussian codebooks, i.e., is drawn according to 
A/'(0,PsJ, for i = 1,2. 

The weighted sum rate maximization (WSRM) problem 
is commonly used to formulate the problem of finding the 
optimal rate region of a multiuser network. However, for the 
two-way AF relay channel, we find two major disadvantages of 
using this formulation. As observed in fi3\, to solve a WSRM 
problem with a specific pair of weights, a bisection algorithm 
is used where each step involves a semi-definite programming 
(SDP) causing high computational complexity. Furthermore, 
from the perspective of practical application, to perform AF 
relaying through this method, every solution corresponding 
to each pair of weights should be obtained and stored in a 
look-up table in the control center before the operation period 
|fT9| . During normal operations, the control center should look 
up the table to decide the appropriate scheme that achieves 
the desirable rate pair Since the channel coefficients may 
vary frequently in different operation periods, the whole look- 
up table should always be updated, which violates real-time 
communication and consumes much more energy. 

To circumvent these disadvantages, we consider an alterna- 
tive formulation of the problem. The basic idea arises from 
the fact that to find an approximate scheme that meets the 
desirable rate pair, we can always regard one transmission rate 
requirement as a constraint. In addition to the power constraint. 



the problem is to maximize the transmission rate of the other 
one under both the constraints. 

Without loss of generality, let's consider the problem of 
maximizing the transmission rate from 5*2 to 5*1 under certain 
power constraint and the constraint of the transmission rate 
from 5*1 to 82- As in the one-way case, we equivalently 
maximize the received SNR at 5*1. The rate constraint is 
converted into the corresponding SNR constraint at S2 as well. 

We consider a general case with both sum and individual 
power constraints, i.e., (Ps^F^ + Ps2^2 + K l) a < 

Psum and al{al + fl^Ps,+ fl^Ps.) < Pk for k ^ 
1, • • • ,71. Then the problem is formulated as 

a^GifafjGia 

a^G2fifrG2a 

a^G2KG2Q; + cr^_2 , (11) 

{Ps.Fl +Ps,Fl+KQi}a< Psu,n 

«^ (^fe + flkPs, + flkPs,) Ika < Pk 
fc = 1, • ■ • ,n 

where Gi — diagjgi}, F; = diagjf;} for i — 1,2, and 71 
is the equivalent SNR constraint of the rate requirement from 
S'l to 5*2. 

Then by Proposition 1 , it is equivalent to solve 
max fjGi/3 

f3,v 

s.t. /3'^GiKGi/3 + erg - 1 < 
/3^Ai/3 - /3^bibf /3 + civ^ < 

{Ps.Fj +Ps,Fl+KQl)f3- PsurnV^ < 

^'^ {^l + flkPs, + flkPs,) Ife/3 - Pkv^ < 
fc = 1, • ■ • ,n 

(12) 

where Ai = 71G2KG2, bi = ^/P^G2iu ci = jia^^2- 

Since mi = 1, by Proposition 2, it follows that 0^ Ai(3 — 
0^hi\y[j3 + c-iv^ < is separated into two convex second- 
order cone constraints given by (/3-^Ai/3 + ciu^) ^ < b^/3 
and (/3^Ai/3 + ciu2)^ < -hi (3. Therefore, two SOCP 
problems are used to solve (fTZt . 

On one hand if the resultant maximum rate satisfy the design 
requirement, then the relays adopt the corresponding optimal 
solution during the data exchange. So the scheme is computed 
instantaneously rather than searching a look-up table. On the 
other hand, the optimal rate region can also be characterized 
via such an approach. This goal can be achieved as follows. 
For the one-way relay channel, the optimal transmission rate 
can be found and it provides an upper bound on individual 
rate for the two-way relay channel. It is equivalent to greedily 
maximizing one transmission rate disregarding the other For 
each rate requirement of one user no larger than the value, 
we may obtain a rate pair which is a boundary point of the 
union rate region. Then the optimal rate region can be fully 
characterized after taking the closure of the convex hull of all 
the derived rate pairs. 



V. Conclusion 

In this paper, we studied both one and two-way AF relay 
channels. Based on a unified formulation, the optimal AF 
schemes for both channels are obtained via numerical algo- 
rithms efficiently. Our future work is to generalize the existing 
results to the complex channel scenario. 
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Appendix A 
Proof of Proposition 1 

Proof: Consider first the following problem. 

x'^bob^x 



s.t. X Aqx + Co < 



(P3) 



c^A,x - x^bjbfx - 



x^AjX + c,j < 0, j 



c,; < 0, 



mi 



1, 



1, • • -TTll 

, mi + 7712 



We show that (PI) and (P3) are equivalent. It is observed 
that the equality should be achieved for the first constraint of 
(P3) at the optima. This can be proved by contradiction. Let 
(xq, Wo) be an optimal solution of (P3). Suppose x^AoXo + 
Co < Wq, it can be found another feasible solution {xo,w'q) 



with Xq AoXo 



Co 



Wn 



such that the objective function 



value of (xo, Wq) is strictly larger than that of Xq AoXo + co < 
Wq. It gives a contradiction of the optimality. 

Then we prove that the two problems have the same optimal 
value. On one hand, let xq denote the optimal solution of (PI). 
It is easy to verify that (xo,w;o) where Xq Aoxo + cq = Wq, 
is a feasible solution of (P3). Then we show that {xo,wq) 
is optimal for (P3). Otherwise, let {x.q,Wq) be an optimal 
solution of (P3), then we get 

xjbobjxo x.'^hoh'^x.Q 



Xn AoXo 



Co 



> 



(xf.)^bob^x[, 
(x[,)^Aox^ +co 



(13) 



where the third equality follows from the previous observation, 
and the inequality and the first equality follows from the 
assumption. Thus it yields a contradiction of the optimality 

of (xf,, 

On the other hand, let (xo, wq) denote the optimal solution 
of (P3). It is clear that xo satisfies all the constraints in (PI). 
Then we show that xq is optimal for (PI). Otherwise, let Xq be 
an optimal solution of (PI) and take = (xq)-^AoXq +co, 
then we get 

x^'bob^xo _ xjbobj^'xo 



xj5 AqXo + Co 



> 



(x[,)^bob^x[, 



^1 \T , 



(14) 



w'q^ (x[,)rAox[, + co' 

where the first equality follows from the previous observation, 
and the inequality and the second equality follows from the 
assumption. Thus it yields a contradiction of the optimality of 

Then we use the trick of variable changing to recast (P3). 
Let y = x/ui and z; = we get 

max y^bobj^y 



S.t. Aoy + cqv' - 1 < 



y' A,y - y' b,b,/ y + cv' < 0, i 
mi + 1, • 



y^A^y + cjv"^ < 0, j 



1, • • - mi 

, mi + m2 

(P4) 



Finally, we prove that to solve (P4), it is sufficient to solve 
(P2). Denote the optimal solution of (P2) as {yopt,Vopt)- 
It is easy to verify that b^yopf > because (0, 0) is a 
feasible solution. Suppose it is not the optimal solution of 
(P4). The rest thing is to show a contradiction. Denote the 
optimal solution of (P4) as {y'opt:'^'opt) ■ Note that (P2) and 
(P4) have the same feasible region. As a result, {y'optiv'opt) 
is also feasible for (P2). By the assumption, it follows that 
i^oyopt) < (boYopt) ' which impHes either h'^yopt < 
^oYopt or ^oVopt < -'^oYopf If the first case holds, then we 
have obtained the contradiction. Otherwise it is easy to check 
that (-y^ 



opt ' opt 



is also an optimal solution of (P4) and is 
feasible for (P2). Consequently, it also yields a contradiction 
to the optimality of {yopt,Vopt)- 

Then we complete the proof. ■ 



Appendix B 
Proof of Proposition 2 

Proof: It is sufficient to consider v > when solving 
(P2). The reason is that only the absolute value of v affects 
the feasibility of the solutions. 

It is observed that the first constraint of (P2) can be rewritten 

as 



[ 



Qo 



y 

V 



1 < 0, 



(15) 



where Qo = diag {Ao, co}. It is a convex quadratic constraint. 
The second kind of constraint can be regarded as the combi- 
nation of two parts shown as follows. 



IQ? 
IQ? 



y 

V 

y 

V 



|2 < hfy, 



< -bf y, 



(16a) 
(16b) 




where Q,; = diag {A,;, c^} is positive semidefinite for i = 
1, • • • , mi. Both of them are convex. Finally consider the third 
kind of constraint. They can be rewritten as 



for Cj < 0, (17a) 
< 0, for Cj > 0, (17b) 



where Qj = diag{Aj,Cj} is positive semidefinite for j = 
mi,-- - ,mi + m2- Note that (17a) is a second-order cone 
constraint and (17b) is a convex quadratic constraint. 

Since cj is a constant, either (17a) or (17b) is considered. 
However, both (16a) and (16b) should be considered. There- 
fore, to solve (P2), one need to solve totally 2"'^ subproblems 
with each constraint in the second kind either recast as (16a) 
or (16b). Note that each subproblem is an SOCP problem and 
thus can be solved in polynomial time. It follows that (P2) 
can be solved in polynomial time because mi is a constant 
independent of the number of variables. ■ 



Appendix C 
Proof of Proposition 3 

Proof: Consider the case when the relay noise dominates 
the sum noise at the destination node. We get 

ctf H2KH2Q;i + 0-2 



R{oLi)= 0.25 log 1 



> 0.25 log 1 + 0.5 ^ ^ ^ 



> 0.25 log 1 



(b) 



afHsKHaai 
ccfHahihfHaaiPs 



afHsKHaai 



- 0.25 



>i?(aopt)-0.25, 



(18) 



where (a) follows from the definition that ccf H2KH2Q;i > 
(T^, and (b) follows from the fact that cxi satisfies the individ- 
ual power constraint and maximizes 

a^HzhihfHzaPs 



SNR^^i (a) 



ct^HaKHaa 



It is easy to verify that SNR^ipi {oti) provides an upper bound 
to SNR{aapi). 

Next, we consider the case when the destination noise 
dominates the sum noise at the destination node. 

R 0:2 = 0.25 log 1 + y —7 

V H2KH2a2 + o-^ / 

> 0.25 log 1 + 0.5^ — I ^ 



( a|^H2hihfH2a2P5 
25 log 1 + — i 



> 0. 



(6) 

>P(aopt)-0.25, 



0.25 



(19) 



where (a) follows from the definition that q;|^H2KH2Q;2 < 
cr^, and (b) follows from the fact that 0:2 satisfies the individ- 
ual power constraint and maximizes 

SNRup2 (a) = 5 . 

It is easy to verify that SNRup2 ((12) provides an upper bound 
to SNR{aopi)- 

Then we complete the proof. ■ 

Appendix D 
Proof of Proposition 4 

Proof: For fairness in comparison, we assume that the 
marginal distribution of the relay noise is the same for all the 
scenarios. Thus if K is a general noise covariance matrix, then 
K I is the noise covariance when noise is uncorrelated. To 
prove the proposition, we compare the following scenarios. 

1) Scenario a: The relay noise is correlated and relays are 
aware of the correlation. Therefore, the relays adopt 
the optimal AF scheme obtained above. Denote the 
corresponding SNR value by SNRa = SNR{cXopi)- 

2) Scenario b: The noise is uncorrelated. Therefore, the 
relays adopt the optimal AF scheme obtained in ifTTl 



denoted as ckq, which is given in the sequel. The 
corresponding SNR value is denoted by SNRi, = 
SNRq (ao), where 



SNRo (a) 



a^H2hihfU2aPs 
aT (H2KH2 I) a + (t2 ■ 



3) Scenario c: The noise is correlated and relays are 
unaware of the correlation. In such a scenario, the relays 
adopt ao- The corresponding SNR value is denoted by 

SNR^ = SNR{cxo). 

If the performance of aopt under correlated relay noise is 
better than the performance of cio under independent relay 
noise, then we say that correlation helps if the relays are 
aware of it. Similarly, if ccq has better performance under 
independent relay noise than under correlated relay noise, then 
we can say that correlation hurts when the relays are unaware 
of it. 

We first revisit ccq as follows. As observed in fTTI . after 
finding the optimal active sets A4*, the optimal AF scheme is 



ao = c (H2KH2 I + (T^G) H2hi , 



(20) 



where G = diag {gi, • • ■ , gn}, gj = 0, j ^ M*, gi, i G M* 
and c are given as follows. 



ieM 



i.mayi i.D i 



at + (J 



E Oti,niax\hi.Dhs,i\ 



\K,Dhs,i\ 



i e M* 



The corresponding SNR value of ao under independent noise 
is given as 

SNRb ^ SNRo (ao) 

= hfH2(H2KH2 0l + a2G)"'H2hiPs, (21) 

It should be noticed that the choice of the optimal active 
set A4 * depends on the network parameters including channel 
coefficients hi, h2 and transmit power constraints Pg and 
Pfe, for k = 1, • • ■ ,n. Nevertheless, it can be easily shown 
that given the amplitudes of all the channel coefficients, Ai* 
remains constant while changing the signs of them. It is 
because in the absence of correlation, K is replace by K I 
and thus diagonal. So, the sign of the amplification gain at 
the relay node k is always chosen such that it coincides 
with sign {hs,khk,D) and the noise power remains unchanged. 
Therefore, the sign of channel coefficient does not affect the 
choice of the optimal active set. Then we get 



SNRc 



a^H2hihf H2aoP5 
a^H2KH2ao +cr2 

aSNRb 



(22) 



where a = a.Q (H2KH2 I) cio + cr^ and 



b = al (H2KH2 - H2KH2 I) ao 



N N 

4=1 j^i 



hl^K{t,j)hl^ 



(23) 



It is easy to check that the second derivative of SNRc 



with respect to b is 



2aSNRh 
(a+b)'-' 



which is always positive since a, 



{a + b) and SNRh are always positive. Thus, SNRc is convex 
with respect to b. Since we have assumed that all the channel 
coefficients are i.i.d Gaussian, sign (/ig.;) and sign{hs,j) are 
independent and take values 1 and —1 with equal probability. 
Then from (|23l l, it follows that the expected value of b over 
all signs of the channel coefficients is conditioned on the 
absolute values of the channel coefficients, while the values 
of a and SNRi, remains constant. Now consider the expected 
value of SNRr 



E [SNRc] = EE.ig 



aSNRb 



(a) 
> E 



aSNRb 



Sn [b] 



E [SNRb 



(24) 



where (a) follows from Jensen's inequality. 

Let us now consider Rc ^ 0.25 log (l + "'^^f" ), the 
achievable rate of otQ under correlated relay noise. It can 
be shown that the second derivative of Rc with respect to 



1 aSNRb{2{a+b)+aSNRt) 
4 In 2 (a+b+aSNRtf {a+b)'-' 



^ is 4i„2 U 

Therefore, Rc 

E[R, 



which is always positive, 
is convex in b. Now consider the average rate 



E[Rc 



EE, 



(a) 

> E 



0.25 log 
E [0.25 log (1 



0.25 log 1 



-iSNRb\ 
a + b J 
^ aSNRb 

a + Esign [b 
SNRb)] = E [Rb 



(25) 



Then we complete the proof. 



